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Subscripts 


ent  = entry 

0 = zeroth  order 

1 = 1st  order 

e = edge  of  the  boundary  layer 


n = surface  normal  direction 
m = represents  mass  transfer 
i,j,k  = characterize  the  finite-difference 
grid  at  the  point  y,  z 


* = iteration 


Superscript 

~ = normalized  value 


INTRODUCTION 


Laminar  incompressible  flows  in  straight  two-dimensional,  axisymmetric 
and  rectangular  channels  have  been  previously  investigated  by  a variety  of 
analytic  and  numerical  techniques.  These  are  typified  by  the  linear  boundary 
layer  (Oseen)  approximation  [lJ  for  evaluating  the  axial  velocity  and  pressure 
history  downstream  of  an  initial  entry  region,  more  exact  numerical  analyses 
using  the  boundary  layer  [2,  3]  or  Navier-Stokes  equations  [4-6]  , and  finally 
a boundary  layer /potential  core  expansion  method  L7-9]  that  models  the  flow 
in  the  entry  region. 

From  these  investigations,  the  internal  flow  in  a three-dimensional  rec- 
tangular channel  can  be  portrayed  by  the  various  regions  depicted  in  Fig.  1 . 


Fig.  1 Channel  geometry 

At  the  leading  edge,  a small  Navier-Stokes  entrance  region,  where  x=0(R  "*), 
develops  R = U—a/Vis  the  Reynolds  number  based  on  the  channel  half-  e 
height  a,  the?  kinematic  viscosity  v , and  the  undisturbed  stream  U . The 
axial  distance  x is  normalized  with  a . Complete  Navier-Stokes  solutions 
are  required  in  order  to  describe  this  entrance  flow.  However,  for  large 
values  of  , this  length  is  so  small  that  extraordinary  axial  grid  resolution 
is  required  ; moreover,  the  surface  viscous  layers  are  so  thin  that  an  ex- 
tremely fine  mesh  is  also  necessary  normal  to  the  channel  walls.  For  these 
reasons,  an  accurate  description  of  this  leading  edge  flow  is  quite  difficult  to 
obtain.  Satisfactory  solutions  are  not  currently  available,  and  this  flow  will 
not  be  considered  here. 

Following  the  Navier-Stokes  region  is  the  so-called  entry  flow,  where 
x - 0(1).  For  this  region,  the  axial  gradients  in  the  surface  viscous  layers 
have  diminished,  and  therefore  boundary  la/er  theory  is  now  applicable.  As 


the  boundary  layers  grow  downstream,  the  potential  core  will  be  appreciably 
affected  by  the  viscous  displacement  effects.  A systematic  matched  viscous/ 
inviscid  expansion  procedure  is  required.  This  analysis  has  been  completed 
for  two-dimensional  cascades  [7]  and  inlets  [8 , 9 J . Agreement  with  solutions 
of  the  complete  Navier-Stokes  equations  [5]  is  quite  good,  so  that  the  validity 
of  the  boundary  layer  procedure  has  been  established.  In  the  downstream 
asymptotic  portion  of  the  entry  flow,  where  x »1,  the  solutions  assume  a 
particularly  simple  form  in  terms  of  the  coordinate  - ^(2x/R  )*.  The  con- 
stant B = 1.  21678,  and  % is  small  but  finite;  e.  g.  , ?<0.  2 . Therefore,  for 
large  R the  entry  flow  can  extend  many  channel  widths  downstream.  In  this 
asymptotic  entry  region  the  axial  velocity  is  uniform  across  the  potential  core 
and  the  effects  of  different  initial  conditions  are  no  longer  discernible  [7-9J  . 

Further  downstream  the  surface  boundary  layers  become  so  thick  that  they 
fill  a significant  portion  of  the  channel.  The  boundary  layer /inviscid  core  ex- 
pansion procedure  is  no  longer  convergent  and  a fully  viscous  channel  flow 
must  now  be  considered.  In  this  flow  domain  ? = 0(1)  and  the  axial  flow  gra- 
dients are  small,  but  finite,  throughout  the  channel.  Boundary  layer-like 
equations  govern  this  flow  regime.  Numerical  solutions  are  required.  Finally, 
for  % - % >>  1,  the  axial  gradients  vanish  entirely  and  fully  developed  con- 
ditions are  achieved.  The  distance  % . has  been  termed  the  channel  "entrance 

* . ..  ent 

length”. 

In  the  present  analysis  all  flow  regimes  downstream  of  the  leading  edge 
Navier-Stokes  flow  are  considered  for  rectangular  inlets.  The  paper  is  given 
in  two  parts.  In  Part  I the  entry  region  is  considered  with  a matched  boundary 
layer/potential  core  procedure  similar  to  that  used  in  two  dimensions.  This 
analysis  applies  for  a three-dimensional  cascade  or  what  is  more  appropriately 
termed  a baffle  configuration.  The  solutions  for  the  asymptotic  entry  flow, 
where  the  effects  of  the  external  geometry  vanish,  [9]  are  obtained  by  trans- 
form methods.  As  in  two  dimensions  a uniform  axial  flow  is  found  in  the  po- 
tential core  region.  Second-order  boundary  layer  solutions,  with  and  without 
surface  mass  transfer,  are  obtained  for  both  the  axial  and  secondary  flows. 

A simplified  result  is  found  for  the  secondary  flow  behavior  in  the  potential 
core;  in  the  boundary  layers,  the  secondary  flow  closely  resembles  that  found 
along  an  infinite  axial  corner  intersection.  [lO.ll] 

In  Part  II,  numerical  solutions  are  obtained  for  a square  inlet,  although  f 
the  analysis  is  developed  for  a rectangular  channel.  These  results  can  be  con- 
sidered to  apply  from  the  asymptotic  entry  region  to  the  fully  developed  region. 
Comparisons  with  the  analytic  solutions  for  the  asymptotic  entry  flow,  where 
the  analyses  of  Parts  I and  II  overlap,  are  possible.  Of  particular  interest  in 
this  investigation  is  the  evaluation  of  the  axial  shear  stress  and  the  accurate 
resolution  of  the  secondary  flow  behavior.  Previous  numerical  studies  [2-4] 
of  this  geometry  have  concentrated  on  the  axial  velocity  and  pressure  distri- 
butions. One  of  the  conclusions  of  the  present  investigation  is  that  this  axial 
flow  behavior  is  virtually  insensitive  to  the  treatment  of  the  secondary  flow, 
although  the  inverse  is  not  true.  In  the  present  study,  converged  secondary 
flow  solutions  are  obtained,  and  uithin  the  limitation  of  the  relatively  course 
grids  considered  here  the  solutions  appear  to  be  quite  reasonable. 

It  is  shown  that  a system  of  boundary  layer  equations  describe  the  entire 
channel  flow  downstream  of  the  entry  region,  where  0(1).  With  an  appro- 
priate transformation  of  the  axial  coordinate  from  x to  ? the  Reynolds  number 
R can  be  scaled  out  of  the  equations.  Reynolds  number  dependence  appears 
oifty  when  the  physical  axial  location  x iT  specified.  A ~ system  of  equations 
including  the  axial  velocity,  axial  vorticity  and  a pair  of  stream  functions  for 
the  rotational  secondary  motion  is  developed.  The  axial  pressure  variation 
or  surface  shear  stress  is  fixed  by  mass  conservation.  A predictor-corrector 
method  is  used  for  the  axial  velocity  and  vorticity,  and  a point  relaxation  pro- 
cedure is  used  for  the  secondary  motion. 

Complete  numerical  solutions  with  and  without  secondary  motion  have  been 
obtained.  Near  the  entrance  a small  overlap  region  where  both  the  analytic 
entry  solution  and  downstream  numerical  solution  apply  is  discerned.  The 
agreement  is  quite  reasonable,  even  though  only  a few  grid  points  lie  within 
the  very  thin  boundary  layers.  The  numerical  solutions  for  the  secondary 
motion  exhibit  a reversed  profile  predicted  by  the  entry  solutions  of  Part  I. 
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PART  I - ENTRY  ANALYSIS 
S.  G.  Rubin  and  S.  Saari 


1. 1 Flow. Description 

The  flow  behavior  in  the  entry  region  of  Fig.  1,  where  x = 0(1),  is  con- 
sidered. The  entry  flow  is  characterized  by  the  formation  of  thin  surface 
viscous  layers  with  a large  central  inviscid  core.  A viscous /iaviscid  ex- 
pansion procedure  is  applicable.  For  the  entry  analysis,  rectangular  channels 
formed  by  a doubly  infinite  baffle  configuration  is  assumed.  A uniform  ir- 
rotational  upstream  flow  is  also  specified.  In  reference  [12]  rotational  entry 
is  also  considered.  The  symmetry  inherent  in  this  baffle  geometry  greatly 
simplifies  the  evaluation  of  the  displacement  induced  core  flow.  The  entry 
analysis  for  a single  inlet  is  very  complex.  For  two  dimensions  this  geom- 
etry has  been  evaluated  by  a Weiner-Hopf  technique  [9];  it  has  been  shown 
that  in  the  asymptotic  entry  flow,  for  x>>1  or  ? = 0(1),  the  effects  of  the  en- 
trance geometry  are  no  longer  important.  Therefore,  the  cascade  and  single 
inlet  solutions  agree  in  this  region.  For  the  rectangular  channels  considered 
here,  the  asymptotic  entry  behavior  in  the  core  region  is  determined  by  a 
Fourier  transform  solution  for  x > >1.  The  entry  flow  for  x = 0(1)  is  de- 
termined by  the  method  of  images.  In  both  cases  the  three-dimensional  rec- 
tangular channel  solution  can  be  represented  by  a pair  of  locally  two-dimen- 
sional solutions  previously  determined  in  reference  [7]  . 

I.  2 Yiscous/Inviscid  Expansion  Analysis- Potential  Core  Solution 

In  the  entry  region  the  large  axial  gradients,  important  in  the  leading  edge 
Navier-Stokes  flow,  are  small  when  compared  with  the  surface  normal  vari- 
ations and  therefore  boundary  layer  theory  is  applicable.  This  behavior  has 
previously  been  observed  in  the  entry  region  of  two-dimensional  channels, 

[7-9]  ; further  confirmation  for  the  asymptotic  entry  region  of  two-dimensional 
inlets,  where  higher-order  terms  in  the  boundary  layer  expansion  are  re- 
quired, is  given  in  section  I.  3 of  the  present  paper.  In  this  section,  the 
leading  terms  describing  the  boundary  layer  and  potential  core  behavior  are 
obtained  for  the  entry  and  asymptotic  entry  flows. 

Boundary  Layer  Analysis 

Thin  boundary  layers 
grow  along  the  walls 
and  an  inviscid  core 
region  is  located  cen- 
trally in  the  channel. 

At  the  corners  the  in- 
terferring  boundary 
layers  lead  to  a 
boundary  region  termed 
the  corner  layer.  This 
flow  has  been  studied 
extensively  for  an  in- 
finite corner,  see  ref- 
erences [10,11],  and 
can  be  neglected  to 
lowest-order  in  the 
boundary  layer  analysis. 
We  will  return  to  the 
corner  flow  later  in 
this  section. 


The  flow  in  the  entry  region  is  depicted  in  Fig.  2. 
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Since  the  flow  gradients  are  large  only  in  the  surface  normal  direction,  the 
standard  boundary  layer  equations  are  obtained  on  each  wall.  JFor  the  re- 
mainder of  the  boundary  layer  analysis  in  this  paper,  the  wall  y = -a  will  be 
considered.  Similar  results  are  obtained  for  the  other  walls.  The  boundary 
layer  equations  are 


u + v + cw  =0 
x y z 

uu  + vu  + Cwu  = - p + — (u  + u ) 

x y z x Re  yy  tt 


p * 0 
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uw  + vw  + oww 
x y z 


-ap. 


+ rr (* 
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+ W ) 

yy 


(I.  la) 
(I.  lb) 
(I.lc) 
(I.  Id) 


where  u,  v,  wj  are  normalized  with  ; x,  y with  a;z  with  b and  p with  eU^  ; 

0 = a/b  = X . The  diffusion  terms  in  the  axial  direction  are  neglected  as 
they  are  0(R  ) smaller  than  the  surface  normal  (y)  variations.  In  the 

boundary  layer,  cross  derivatives  (z)  are  also  small  when  compared  with 
normal  changes,  but  are  important  in  the  corner  layers  and  may  be  required 
for  the  second-order  boundary  layer  analysis  of  the  following  section.  Using 
standard  boundary  layer  assumptions,  we  have 
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(I.  2a) 
(I.  2b) 


w = Rj  ^ Wj  + Rj 1 w2  + . . . , 

where  y = 0(R  ^ ) ; x,  z = 0(1).  To  lowest  order  in  R ^ 
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equations  (1. 1) 

(I.  3a) 
(I.  3b) 
(I.  3c) 


The  axial  pressure  gradient  is  zero,  since  the  inviscid  stream  is  undisturbed 
and  uniform.  Equations  (I.  3)  with  the  boundary  conditions  u = v = o at  y=-l, 
and  u->UB  for  y-ye  are  precisely  those  of  the  flat  plate  boundary  layer.  The 
subscript  e denotes  the  edge  of  the  boundary  layer.  The  solution  of  these 
equations  is  well-known  C 13 J and  given  by 

Uo  " 'o'  (T1)  • *1  = (2xRe>‘i  <*)  fD'  <n)  * f,0l»  . (I.  4a) 


where  *• 

'o  (T’)  + 'o  in)  C ™ = 0 • 'o  (o)  = fo  (o)  = 0 • ('•  4b) 

( (o)  = 0.  469600 

T)*  yRe^(2x)'i 
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(I.  4c) 


*4 


Introducing  the  notation  ?2  = f?  2x/R  , we  have  T)  = . For  y-y 

tl-»  • . we  find  ^ - 


v P (2xR„) 


p2  ?_1  R "l  , where  0 = 1.  21678  . 
e 


or 


Potential  Core:  Displacement  Effects 

The  surface  boundary  layers  perturb  the  central  inviscid  core.  The 
effective  core  area  is  diminished  thereby  accelerating  the  flow.  Since  the  in- 
viscid flow  is  irrotational  the  JWree-dimensional  Laplace  equation  must  be 
solved  for  the  perturbed  core  fl&.v.  Capital  letters  are  used  here  to  distinguish 
the  potential  flow  variables  from  those  of  the  boundary  layer.  A unified  no- 
tation is  used  in  the  figures.  Therefore, 

V = Um  [x  + ^ ^ (x,  y,  z)J  , and  7^  rp^  = 0 . 

Also, 

U = cp  =U_+R'^U,  . etc. 
x • e 1 

By  matching  with  the  inner  boundary  layers,  we  obtain  for  x>o  the  boundary 
conditions 

T,  =P(2x)'*  (I.  5a) 

n 

on  all  walls,  where  n denotes  the  surface  normal  direction  measured  positive 
inward;  e.g.  on  y=l,  n=  -y,  etc.  For  x<  o,  with  the  baffle  configuration,  the 
symmetry  conditions  lead  to 


= 0 


(I.  5b) 


Finally,  for  x-.-  m we  require  the  uniform  flow  condition,  ^ 
In  view  of  the  symmetry,  and  the  fact  that  the  boundary  * 
are  independent  of  y and  z,  we  can  write 

vpj  (x.  y.  z)  = V j (X,  y ) + tp2  (x,  z)  , 

1 


= 0. 

conditions  (I.  5) 


where  <p.  and 
conditions  on 


satisfy  two-dimensional  Laplace  equations.  The  boundary 

y-  ±1  apply  to  cpj  , since  :pj  • 0,  and  those  on  z = *1  apply 

to  *p^  . Therefore  the  three-dimensional  Laplace  equation  and  associated 

boundary  conditions  are  satisfied  by  a pair  of  two-dimensional  solutions  each 
corresponding  to  the  entry  solution  in  a two-dimensional  cascade.  This  is  a 
significant  simplification,  since  the  two-dimensional  solutions  are  given  in 
Van  Dyke  [7]  by  the  method  of  images  and  by  Fourier  transform  techniques. 
The  complete  three-dimensional  analysis  and  details  of  the  solution  are  given 
in  reference  [12].  The  image  technique  is  preferable  for  the  entry  flow,  while 
the  transform  procedure  is  preferable  in  the  asymptotic  entry  region.  The 
analysis  and  solutions  are  summarized  below; 

Image  Solution 

With  the  appropriate  image  or  Green's  function  distribut  ons,  we  obtain 
the  following  solutions: 

v’,  • i ifr-VY /•  i ( v?:: )* 
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Note^that  Cp^  does  not  appear  in  this  expression. 

W.  = cp  is  obtained  from  (I.  6a)  with  y — Oz  and  n^*an  inside  the 
summation  n = 1 to  •. 

* jpw-* 

+ f<-4-+  <2n-y )f  )*-»y  + /(«f,-K2nty)V-x\  ‘1 

\ * + (2n-y)  / \x2  + (2n+y)2  / J 


-(J^4!  * ip" 

, Ax2  + (2m- z)2  02^-xV  Ax2  + 

\7T7ii^ri?  ; \ir 7 


(m-l)‘  -m5 


2 2 ^ 
0*  (Zmlzf)  • 

a2  (2m -fz)^ 


The  centerline  axial  velocity  for  different  rectangular  inlets  is  shown  in  Fig.  3 

for  R = 75.  Typical  solutions 

— , T . , for  the  streamwise  isovels  at 

“t  At . 79  various  locations  in  the  entry 

. > flow  are  shown  on  Fig.  4 

(J.  for  different  rectangular  tn- 

lets.  It  is  interesting  to  note 

that  for  very  small  values  of 

is-  ^t^**0**’  x behavior  is  monotonic. 

F°r  larger  x values  a valley 
\00‘^0~*00‘  *n  axial  flow  distribution 

i4-  y is  Predicted*  This  behavior 

^ I has  been  observed  in  two- 

dimensional  Navie  r- Stokes 
solutions  [.4-6].  Finally,  for 
even  larger  x values  the  be  - 
havior  is  once  again  mono- 

1^—T\ j a.  * . tonic,  but  the  direction  of  in- 

**  0 1 * * 4 • ■ creasing  u has  changed  from 

°*,K>  ow*  i that  found  for  small  x.  A 

typical  solution  fftr  a 5x1 

Fig.  3 Axial  v#locity  distribution  in  chann  1 is  shown  on  Fig.  4d. 

entry  region  -rectangular  channel 


Transform  Solutions 

U the  potential  problem  is  solved  by  Fourier  transforms,  the  same  pair- 
wise two-dimensional  behavior  occurs  and  the  solution  for  V.  is  as  follows 

[7.12]  : 1 

m 

ip  (x.y,«).-I»  J l1-*  ■*?!"«)«  ~t"’<  co»h  ■»-/)  d.  (1.7) 

m - *1 »| ■ »inh  * 


[ (1  + i.£n  »)e'i,x  cosh  o.(l-z) 

J .1.7* finlT  n. d* 


Fig  4a  Streamwise  isovels  at  inlet - 
square  channel 


Fig.  4b  Streamwise  i9ovels  in  entry 
region- square  channel 


*••1000  x.io  1*1  o 

(•0  05441 


Fig.  4c  Streamwise  isovels  in  entry 
region- square  channel 


Fig.  4d  Streamwise  isovels  near 
entry-  rectangular  channel 


The  velocities  etc.  are  difficult  to  evaluate  for  general  values  of  x,  but 
for  x>>1  (the  asymptotic  entry  flow)  it  is  possible  to  obtain  asymptotic 
expressions  using  the  generalized  Fourier  transform  concept  outlined  by 
Lighthill  C 14 J . More  complete  details  of  this  procedure  can  be  found  in  ref- 
erences [7,12]  . The  results.for  x>>  1,  lead  to  the  following  expansions  of 
solution  (I.  7). 
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Uj  = ^ = 2*  B Lx^  (1  + 0)  + 


1 -a/a  2 a/>  2 "I 

^ x ' Oy*2  -1)  + C r(3z‘i  -1)  J + . 

Vj  = ^ = -2*^  By  [x~^  + (l-yZ)/8x  ^ J + . 

y 

Wj  = = -2  * Bz[_x'^  + (1-z2)/8x  ^ J + 

z 


(I.  8a) 


The  non-monotonic  behavior  seen  earlier  in  the  series  solutions  is  still 
apparent  in  the  asymptotic  solution.  In  reference  iiZi  the  asymptotic  and 
series  behavior  for  x<o  is  presented.  Here  we  are  concerned  only  with  the 
internal  flow. 

For  very  large  x>>1,  the  full  solutions  for  the  velocities  become 


•U  -1  + (1  + o)?  + 


-1 


-1 


y+ 


W.-bV'r  -1 


(I.  8b) 


Therefore  the  asymptotic  entry  flow  and  onset  of  the  fully  viscous  channel  flow 
have  the  following  properties  in  the  central  core  region: 

(1)  The  axial  velocity  is  uniform  across  the  channel. 

(2)  The  relevant  axial  flow  coordinate  is  % - 0(1)  i 

(3)  The  secondary  velocities  are  0(Re  *)  and  not  0(Re  / as  in  the  entry 
flow  boundary  layers  and  potential  core. 

This  behavior  has  previously  been  found  for  two-dimensional  channels  and 
corresponds  to  the  form  originally  assumed^by  Schlichting  [13]  in  his  expansion 
technique.  For  two-dimensional  channels  V is  given  by  (I.  8b),  W=  0 and 
U~  1 + % . 

Fig.  5 depicts  the  axial  flow  variation  along  the  channel  centerline  as 

obtained  by  the  image 
and  asymptotic  solutions. 
It  is  interesting  that  the 
asymptotic  form  is 
achieved  for  relatively 
small  values  of  x,  even 
though  it  has  been  de- 
rived  for  x»  1.  Similar 
behavior  is  found  at 
other  channel  locations. 
This  allows  for  the  use 
of  the  simplified 
formulas  (I.  8b)  over  a 
significant  portion  of 
the  entry  region. 


Fig.  5 Entry  flow  solutions-  Square  channel 
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Finally,  Fig,  6 depicts  the  qualitative  behavior  of  the  secondary  flow 

stream  lines  projected  onto 
the  (y,  z)  cross  plane.  For 
a square  channel  these  lines 
y are  straight,  for  a 2x1 

channel  parabolic,  and  in 

0| — — 1 general  they  are  given  by 

the  equation: 

dy  ^1  y 1 y 

y//k  inrs^a  * ~ t • 


- - o 

z = cy 


Z : cOy 


where  c is  a different  con- 
stant for  each  of  the  stream- 
lines. For  a square  channel 
z=  cy;  for  0 = 2 , z=2cy^  , 
etc. 


Fig.  6 Asymptotic  flow  projection  on 
cross  plane. 

I.  3.  Second-Order  Boundary  Layers  - Asymptotic  Entry  Region 

With  the  asymptotic  potential  solutions  (I.  8b)  the  second-order  boundary 
layers  can  now  be  evaluated.  With  the  expansions  (1.2)  in  equations  (I.  l),the 
axial  and  secondary  motion  along  the  surface  y = -1  are  described  by  the 

following  equations.  Terms  of  0(R  z ) are  retained. 

e 

u.  + R ^ v,  t 0 w.  = 0 , 

1 e 2y  1 

x 1 z 


u u + u.  u -t  R 


‘ V.  U.  -f  V-  R ’ ' u + p,  = fr—  U, 
e 1 ly  2 e oy  rlx  Rp  1 


u w.  f R 
o 1 e 


*l*ly  " '0pl*  *iT  Syy 


As  y-*y^  , , Wj-*Wj  from  (I.  Hb);  therefore,  -pj  = (l+0)B(2x)  , 

.3/3  .3  /3  x 

and  - Cp.  = €2  B*  x ; € = 1 for  rectangular  channels  and  0 = C = 0 for 

two  -dimensional  channels. 

Assuming  that 
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u = (n)  + (l+a)?fj  (11)  = ao  + Re’i  Uj 

wj=  -£e*(2X)-*h'(Ti)  = Wjh'(Ti) 
the  following  equations  are  obtained  for  fj(T|)  and  h (T) ) : 

h“  + (f  h')'  = 1 . (I.  9a) 

and 

c + f f*,  - f'  l\  + 2f*  f,  = -1  + e f*  h/ (l+o)  , (I.  9b) 

lolol  ol  o 

where  from  continuity  v^  = e j_(  1 + o ) (T)  f ^ - 2fj(r|))  + eh(T))J  • 

The  boundary  conditions  become 

fj  (o)  = f{  (o)  = h'(o)  = 0 , (I.  9c) 

lim  fj  (T))  = h'(T|) -.1 
T)-« 

Equation  (I.9b)  for  the  secondary  velocity  w.  is  identical  with  that  derived  in 
earlier  studies  of  the  flow  along  an  unbounded  corner  [lo3  . An  interesting 
secondary  flow  reversal  is  predicted,  and  this  behavior  is  shown  qualitatively 
on  Fig.  6.  The  core  solution  was  discussed  in  section  1.2  . The  reversed 
flow  profile  is  depicted  more  accurately  in  section  I.  4 . 

The  second-order  axial  flow  is  described  by  equation  (I.  9a).  For  two- 
dimensions  (c  = o)  this  equation  was  considered  by  Schlichting  [13]  and  Wilson 
[8J  . Solutions  are  shown  on  Fig.  7 . Also  presented  is  the  potential  core 
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Fig.  7a  Two-dimensional  entry 
flow  solutions 


Fig.  7b  Two-dimensional  entry 
flow  solutions 


behavior,  the  first-order  Blasius  solution,  and  a composite  viscous /inviscid 
profile.  Comparisons  are  made  with  the  Navier-Stokes  solutions  of  Wang  and 
Longwell  [5]  . The  agreement  is  reasonably  good  for  £ < 0.  2 and  therefore 
the  applicability  of  the  expansion  procedure  in  the  entryTlow  is  confirmed. 

For  three-dimensional  flows,  solutions  have  been  obtained  for  a square  inlet 
(0  = 1).  These  results  will  be  presented  in  Part  II,  where  comparisons  with 
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Fig.  7c  Two-dimensional  entry 
flow  solutions. 


the  viscous  channel  solutions 
are  included.  The  boundary 
layer-like  behavior  found  in 
(I.  8)  is  used  to  describe  the 
downstream  viscous  channel 
flow.  The  expansion  procedure 
shows  evidence  of  divergence 
for  values  of  j >0.2  . Higher- 
order  terms  (not  shown)  exhibit 
behavior  typical  of  divergent 
asymptotic  series.  Also,  eigen- 
values of  equation  (I.  9a)  must 
then  be  included  in  the  expansion. 
These  introduce  unknown  con- 
stants which  must  be  evaluated 
by  other  means,  see  references 

L 8 , 9]  . 

A final  comment  on  the 
secondary  behavior  in  the  entry 
region  concerns  the  corner 
layers  of  Fig.  2 . The  corner 
flow  for  an  infinite  geometry  has 
been  presented  in  great  detail  in 
earlier  studies  [10, 1 1 . Sig- 
nificantly, for  the  channel  the 
secondary  flow  behavior  near 
the  corners  (z  — ±1)  gives 

w.  - ± 9 (2x)  R h'  Cn)  . 

1 e 


This  is  precisely  the  result  found  for  the  unbounded  corner  flow  and  therefore, 
in  the- asymptotic  entry  region,  the  solutions  of  references  L10,  llj  apply  to  the 
channel  as  well.  A swirling  motion  shown  qualitatively  on  Fig.  6 is  obtained. 

This  completes  the  entry  flow  description  for  the  axial  and  secondary 
motion  in  the  potential  core,  the  boundary  layers  and  corner  layers.  Com- 
parisons with  numerical  solutions  are  given  in  Part  II.  In  section  I.  4 a brief 
description  of  the  flow  behavior  allowing  for  surface  mass  transfer  is  pre- 
sented. The  effect  of  injection  is  particularly  interesting. 

1.4  Surface  Mass  Transfer 

In  view  of  the  relative  simplicity  of  the  asymptotic  entry  flow  analysis,  the 
effect  of  surface  mass  transfer,  both  suction  and  injection,  on  the  secondary 
motion  in  this  region  is  now  evaluated.  The  effects  of  mass  transfer  on  the 
streamwise  profiles  have  previously  been  calculated  for  a flat  plate  bv  the 
similarity  analysis  of  Emmons  and  Leigh  Ll5j  . The  mass  injection  is  of  the 
type  v~x  • , or  in  terms  of  % and  the  normalized  velocity  v=*R  9 v , 
we  obtain  v constant.  For  a uniform  outer  stream,  as  occurs  in  the  asymp- 
totic entry  region  for  small  ' , the  modified  boundary  layer  equation  for 

u ('  (ri)  on  the  lower  channel  wall  becomes 
m 

C -t  f i"  = o 
m mm 

whe  re 


f (o)  - m,  f ' (o)  - 


lim 


[f#  (T|  )-.!,  f (T’)~  r\  -r  ] 


r\-n 


The  normalized  surface  velocity  v is  given  as 
v = p ’ m 


m>o  for  injection;  m<o  for  suction,  P = Pq  - 1.21678. 
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The  solutions  for  f are  tabulated  in  reference  L 1 5 J . The  asymptotic 
boundary  condition™  on  the  displacement  velocity  v becomes 


lim  „ v~  0 / 0 

*1-  m m 


9 is  tabulated  for  several  values  of  m 
m 

Table  1:  Shear  stresses  with 

mass  transfer 


m 

(0) 

o°> 

-0.707 

3.8118 

0.0502 

-0.493 

2.3913 

0. 1312 

-0.333 

2.0929 

0.2326 

-0.177 

1 .6108 

0.3460 

-0.0707 

1.3663 

0.4191 

0.0 

1 .21678 

0.4696 

0.0707 

1.0733 

0.5214 

0. 177 

0.8766 

0.6012 

0.333 

0.3724 

0.7394 

0.495 

0.3473 

0 8538 

0 707 

0.0330 

1.0308 

0.778 

- 0.0670 

1.0910 

0848 

-0.1650 

1 . 1518 

in  Table  I;  also  presented  are  values 
of  the  shear  stress  parameter 

fm(o)  * 

The  normalized  secondary 
velocity  w can  be  defined  as 

* ^Bm/s>h'm  (T1) 


where 


h'"  + f h'  + f ' h ' 
m mm  mm 


X 


and 


also 


h ' (o)  = o,  lim  h' 

m _ n 

T|-» 


h'  (o)  = -B 

rr>  1 ' 


(Tl)-l 


The  solutions  for  h#  and  i' 
are  shown  in  Fig.  8.  m 
Mass  injection  increases  the 
secondary  flow  reversal  as  well  as 
the  streamwise  surface  shear 
stress.  With  large  suction  the 
secondary  flow  reversal  can  be 
eliminated.  On  Table  1 it  is  seen 


Fig.  8 Effects  of  mass  injection  on  asymptotic  velocity 
profiles  in  entry  region. 

that  h " (o)  changes  sign  at  about  m=0.71.  The  potential  core  solutions  (1.8) 
for  v and  w are  modified  by  the  factor  p /8»  but  the  projected  streamline 
pattern  is  unchanged.  In  the  surface  bouncfiry  layers  the  secondary  motion  is 
accentuated  or  reduced  as  depicted  in  Fig.  8 . It  is  evident  from  the  figure 
that  the  secondary  flow  patterns  are  significantly  altered  with  even  moderate 
mass  transfer. 
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Part  II.  - Numerical  Solution  For  A Square  Channel 
S.G.  Rubin  and  P.  K.  Khosla 


II.  1 Flow  Analysis  - Fully  Viscous  Region 


The  boundary-layer/potential  flow  expansion  procedure  applicable  in  the 
entry  region  Fig.  1 and  described  in  Part  I fails  in  the  downstream  flow  where 
x>  > 1 or  % = 0(1)  . From  the  results  of  Part  I we  expect  that  the  entry  anal- 
ysis (both  series  and  asymptotic  solutions)  will  apply  many  diameters  down- 
stream, e.  g.  , for  all  § < 0.  20;  but,  for  larger  % values,  the  errors  in  the 
asymptotic  boundary-layer /potential  core  expansion  procedure  increase  as 
the  boundary  layers  grow  and  begin  to  fill  the  channel. 

From  the  results  of  Part  I,  it  is  seen  that  a fully  viscous  model  is  required. 
The  equations  are  appropriately  written  in  terms  o£  the  coordinates  , y,  z). 

The  normalized  velocities  are  u=u/U— , v=?R  v/U^  w=§  Rg  0 v// , 

and  the  normalized  pressure  is  — , 2 , 

p=  p/ 0 u— , where 

p = pQ  + r;1  pi  + r;2  p2 


See  Part  I,  equations  1.8  for  the  asymptotic  forms  required  here.  Rdaining 
only  the  lowest  order  terms  in  ?R  in  each  of  the  governing  equations  (see 
Part  I,  equations  1.1),  we  obtain: 


Continuity: 

? Momentum: 


u.  + v + Ow  = o 

? y i 

58-V.uu,  + vu  + cwu  + p 

% y z ro§ 


(n.i) 

(H.2) 


y Momentum: 

p * O , 

K°y 

Pi  = ° 

Kiy 

(II.  3) 

z Momentum: 

p = o , 

o z 

Plz  = ° 

(II.  4) 

Streamwise  Vorticity:  fi  = 

0V*-Wy 

<n.5) 

?8  ^ = uQ»  ♦ TV,  ■ lW,  • 5 * (ufi  t TV  - aw) 


+ vfl  +o  w 0 + 0 ( v + o w ) , (II.  6) 

y z y z 

2 2 

where  t = ou  , a=  u , y ■ ( ) + 0 ( ) , 

z y ' 'yy  zz 

and  0 = a/b  \ = 1 for  the  square  channel  considered  here.  Equations 
(11.1,6)  are  reformulated  in  terms  of  a potential  functions  and  a stream 
function  ♦ for  the  crossflow  (v,w).  We  define 


v = rp  ♦ c ♦ and 

y * 

w = av  - * 

* y 


(n.  7) 


Therefore,  (II.  1,  5)  become,  respectively, 
2 

y v = -u.  , 

and 

y2  * = o 


(II.  8) 


<n.9) 
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Equations  (II.  2,  6,  8,  9)  describe  the  unknowns  u,fl,cp,  4.  v and  w are  re- 
covered from  (II.  7)  . The  streamwise  pressure  gradient  depends  only  on 
as  seen  from  (II.  3,4)  and  is  specified  b>  the  mass  conservation  integral 

JJ  u (§,  y,  z)  dydz  = 1.  (11.10) 

A 


A denotes  the  rectangular  cross-sectional  area. 

The  appropriate  boundary  conditions  are  shown  on  Fig.  9.  It  is  easily 


proven  that  these  conditions  uniquely  define  cp  and  i)i  . The  zero  normal  mass 
flux  condition  on  v and  w is  satisfied  exactly.  The  no- slip  conditions  are 
satisfied  indirectly  through  the  surface  vorticity  relation. 

For  two-dimensional  channels,  solutions  were  also  obtained  with  an  al- 
ternate formulation  in  terms  of  the  stream  function  and  vorticity.  The  surface 
shear  was  determined  with  the  mass  conservation  condition  (II.  10)  . The 
pressure  gradient  pxwas  then  obtained  by  application  of  the  momentum 
equation  (II.  2)  at  the  surface,  where  the  no- slip  conditions  are  applied.  It  is 
significant  that  fully  converged  solutions  obtained  at  each  streamwise  location 
by  the  two  procedures  were  identical.  The  accuracy  of  the  final  solution  is 
therefore  limited  solely  by  the  mesh  size.  A similar  vorticity  formulation  is 
possible  for  three-dimensional  flows,  although  it  has  not  been  used  here. 

Another  important  feature  of  the  present  formulation  is  the  conservation 
implicit  in  the  system  (II.  8,  9).  With  the  Neumann  boundary  conditions  on  cp 
it  is  necessary  that 


j*  J cp  dydz  = 


(II.  ID 


A 


With  the  mass  conservation  condition  (II.  10)  we  automatically  satisfy 

J j **5  dydz  = ° * (xi.  i2) 

A 

and  therefore  the  integral  condition  (II.  11)  is  always  satisfied.  In  an  earlier 
study  of  this  problem  Briley  [2 J applied  a similar  procedure  with  a single  po- 
tential function,  but  it  was  necessary  to  apply  an  iterative  correction  factor  in 
order  that  (II.  11)  be  satisfied. 

More  important,  however,  the  additional  stream  function  provides  a simple 
mechanism  for  exactly  satisfying  the  vorticity  (0)  and  secondary  flow  (v,w) 
equations  to  any  desired  degree  of  accuracy.  In  this  way  we  are  able  to 


The  zero  subscript  on  p has  been  dropped. 
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I' 1 


rigorously  examine  the  mutual  interaction  between  the  streamwise  and 
secondary  motions  and  thereby  accurately  evaluate  the  effect  of  one  on  the 
other.  Finally,  the  asymmetry  of  0 leads  to  the  additional  condition 

r p 2 

J J V i-  dydz  = o 

A 

This  is  insured  by  the  boundary  conditions  relating  i|i  and  0 . 

The  governing  system  for  the  downstream  flow,  % - 0(1),  and  ultimately 
the  fully  developed  region,  ?>>  1,  has  now  been  described.  It  is  significant 
that  the  Reynolds  number  R does  not  appear  directly,  either  in  the  equations 
or  in  the  boundary  condition^.  The  Reynolds  number  dependence  exists  only 
in  the  transformation  relating  % to  the  physical  distance  x.  Therefore,  with 
the  present  analysis,  the  Reynolds  number  does  not  enter,  in  any  way,  into 
the  numerical  calculations.  A single  solution  valid  for  all  Reynolds  numbers 
is  obtained.  If,  for  example,  we  consider  a channel  defined  by  0 <§<1,  then 
the  Reynolds  number  defines  the  physical  length  of  the  duct.  As  fF  Tncreases 
the  duct  length  increases  and  vice-versa. 

Finally,  an  overlap  region  where  both  the  entry  solution  of  Part  I and  the 
numerical  downstream  results  are  applicable  should  exist.  We  recall  from 
Part  I that,  for  a square  channel,  the  asymptotic  behavior  found  in  the  po- 
tential core  of  the  entry  region  is  of  the  form 

u ~1  + 2?  , v-  -y,  w ~ -z  , (II.  13) 

while  the  boundary  layer  variation,  away  from  the  corner  regions, is  as  follows: 

u ~i'Q  (T))  + (T))  , (11.14) 

w * - zh  ' (n  ) 

These  solutions  approximate  the  viscous  channel  equations  (11.1,2,6)  for 
small  values  of  ? . Therefore,  the  numerical  solutions  for  the  downstream 
viscous  channel  flow,  which  are  valid  only  for  x >>1,  should  become  meaning- 
ful for  moderate  values  of  § , where  (II.  13, 14)  are  indicative  of  the  flow  be- 
havior. The  complete  channel  flow,  away  from  the  small  leading  edge  Navier- 
Stokes  region,  is  then  described  by  the  composite  entry  plus  downstream 
solution. 


II.  2 Numerical  Methods 

The  governing  equations  (II.  2,  6,  8,  9, 10)  are  solved  numerically  by  a 
marching  procedure  starting  at  the  leading  edge  % = o.  The  initial  conditions 
are  assumed  to  be  Cl  =v  =w  = o.  The  axial  velocity  u is  uniform  and  satisfies 
the  mass  conservation  condition  (II.  10).  At  the  walls  the  no- slip  and  zero 
mass  transfer  conditions  are  enforced.  Since  only  the  quarter  channel  is  con- 
sidered, symmetry  conditions  are  prescribed  along  the  center  lines.  The 
complete  set  of  boundary  conditions  is  shown  on  Fig.  9.  These  initial  values 
are  somewhat  artificial  and  therefore  the  numerical  solutions  become  mean- 
ingful only  after  the  influence  of  the  initial  values  become  negligible,  i.e.  , in 
the  asymptotic  entry  flow  described  in  Part  I.  In  view  of  the  parabolic 
character  of  the  viscous  channel  equations,  the  only  remaining  influence  of 
the  initial  values  will  be  reflected  in  the  mass  conservation  integral  (10).  This 
behavior  has  previously  been  verified  for  two-dimensional  channels  [4-6]  , 
where  the  asymptotic  entry  behavior  is  closely  approximated  by  the  numerical 
boundary  layer  solutions  [7-9]  . As  will  be  shown  in  the  present  analysis,  a 
similar  overlap  is  found  for  a three-dimensional  geometry. 

Finite  difference  discretization  is  used  for  all  derivatives.  These  are 
centered  at  the  half  point  ?.  x in  marching  from  %.  to  . ; i.  e.  , a Crank- 
Nicolson  procedure.  Therefore,  1 11 


. “i+i,  j.k‘ui,  i.k 

i+i.i.k TT ' 


(II.  15a) 


(II.  ISb) 


f 


. > = ui,j+l.k'ui,j-l,k  + Ui+1,  j+l,k~Ui+l,j-l,  k 

' y'i+t.j.k  IH 


. . _ Ui,  j+1,  k~2ui,  j,  kiui,  j-1,  k*ui+l,  j+1,  k~2ui+l,  j.  k+ui+l,  j-1,  k 

' yy’i+i,  j.  k ^2 


(1115c) 


All  other  variables  are  discretized  in  a similar  fashion;  z derivatives  are 
obtained  from  (11.15)  with  j-k,  k-j  and  then  transposing  the  subscripts. 
Non-linear  terms  such  as  uu,  in  (II.  2)  are  treated  with  quasi-linearization, 
so  that  * 2 *2 

(uu, ) , . 2Ui+l.).k  Ui+l,j,k  -“i.j.k^i-tl.j.k 


where  6?  denotes  the  step  size  and  the  asterisk  superscript  denotes  the  value, 
at  the  new  location  , , obtained  from  the  previous  iteration.  The  function 
u(?:.yr  \)  i»  denoted  1 by  u where  ? = (i-l)6?,  y«  =(j-l)h  and  z . ={k-l)h; 

h fs  1 the  uniform  grid  spacing  in  the  (y,  zf  plane  so  that  h=l/N.  N is  the 
number  of  grid  intervals  in  the  quarter  channel,  see  Fig.  9. 

For  each  equation  (II.  2,  6,  8,  9)  we  obtain  a system  of  the  form 


A.  u...  ...  ,+B.  u.  . . . , +C.  u. , + D.  u.4l  . ... 
l l+l.j+l,  k l l+l,  j-1,  k i i+I,  j,k  l l+l,  j,  k+1 


(II.  Ik) 


+ E.  u.Al  . . . = F. 
l l+l,  j , k- 1 i 


A F.  depend  only  on  the  values  at  ?.  or  of  the  previous  iteration 

(asterisks1).  The  full  equations  are  easil^r  derived.  Equation  (11.16)  is  solved 
with  an  ADI  predictor-corrector  method  for  u,fi  and  point  SOR  for  Cft  ♦ . For 
the  former  the  tridiagonal  matrices  are  inverted  with  the  now  standard  two 
pass  algorithm.  In  the  ADI  iterative  procedure  the  u.  . . . , , u.  . . . and 
u.+i  . j k terms  are  treated  implicitly  in  the  predictor’  * ’ 1 > J’  stage 

1 ’ 1"  ’ and  the  terms  u. . , , ,,  u.  . , . , and  u.„  . . .are  implicit  in 

the  corrector  stage.  This  l+1’  >’  k+1  l+1>  >’  k 1+1*  k procedure  was 

found  preferable  to  the  line  predictor-corrector  method  [16]  as  the  symmetry 
properties  of  the  flow  variables  were  maintained  to  a high  degree  of  accuracy. 
For  the  SOR  calculations  an  optimum  over- relaxation  parameter  UJ  Cl7 D for  a 
square  uniform  grid  with  Dirichlet  boundary  conditions  is  specified; 
id»2  (1-tt/(N+1).  This  parameter  also  works  quite  well  on  the  Neumann  pro- 
blem (II.  8)  for  . For  N=10  1.  43.  For  N=20  utiwl.  72.  For  all  equations, 

the. system  (11.16)  is  diagonally  dominant  and  appears  to  be  stable  even  for 
relatively  large  values  of  . 

Finally,  the  pressure  gradient  in  (II.  2)  is  updated  by  Newton  extrapolation 
such  that  the  mass  flux  condition  is  satisfied.  The  mass  flux  integral  is 
approximated  with  the  trapezoidal  rule.  The  secondary  vejocities  v and  w 
are  recovered  from  (II.  7). 

The  iterative  procedure  can  be  summarized  as  follows: 

At  station  ?.+j  let  p,  assume  their  ?.  values. 

Calculate  u.+j  from  (II.  2)  using  Iterative  ADI. 

Calculate  (l+j  from  (D.6)  using  interative  ADI. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


and  are  found  from  (11.8,9)  using  SOR. 


*i+l 

Update  Vj+j  and  w.^,  from  (II.  7) 


i+1 


Update  (P^)j^j  from  (II.  10)  and  return  to  step  (1). 


At  each  £ location  this  procedure  continues  until  all  values  have  converged 
to  the  specified  tolerance,  e.g.  six  decimal  place  accuracy. 

In  the  vorticity  streamfunction  method  discussed  earlier,  the  evaluation 
of  p^  in  step  (6)  is  replaced  with  the  direct  evaluation  of  the  wall  shear,  e.g. 
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T— U—  on  y=-l,  such  that  the  man  flux  condition  (II.  10)  ia  satisfied.  The 

firessure  ia  then  obtained  by  application  of  (II.  2)  at  y--l.  The  two  procedures 
ead  to  identical  results  in  two'dimenslonal  channels.  This  latter  approach 
would  be  preferred  for  problems  in  turbomachinery  where  local  separation 
bubbles  can  be  evaluated.  This  approach  has  previously  been  considered  only 
for  external  flows  [18]  . 

II.  3 Results 

The  numerical  solutions  were  obtained  by  considering  the  quarter  channel 
with  mesh  widths  h = 0.05  (21  points)  and  h = 0. 1 (11  points),  respectively. 
Therefore,  in  the  entry  flow  where  the  boundary  layers  are  very  thin,  the 
accuracy  of  the  results  near  the  channel  walls  is  marginal.  Further  down- 
stream, as  the  boundary  layers  thicken,  the  resolution  increases  markedly. 

All  of  the  results  are  jfiven  a,s  functions  of  the  normalized  axial  length  ? , 
where  we  recall  that?”  = 2x(T  / R • The  solutions  are  independent  of  the 
Reynolds  number,  except  for  the8 normalization  inherent  in  ?,  v,  w.  NoleThat 
for  large  R even  relatively  small  values  of  ? can  correspond  to  axial  dis- 
tance of  ma^y  channel  widths. 

The  centerline  velocity  is  depicted  on  Fig.  10a  for  the  entry  flow,  and 
Fig.  10b  shows  the  entire  channel  flow  to  fully  developed  conditions.  Solutions 

with  the  coarse  and  fine 
grids  are  presented,  as 
are  solutions  with  zero 
crossflow  (secondary 
flow).  The  different 
initial  conditions  with 
differing  meshes  result 
from  the  requirement  that 
the  mass  flux  conditions 
(II.  10)  be  satisfied.  Also 
shown  are  the  entry  so- 
lutions obtained  in  Part  I. 
For  a Reyholds  number 
R =75  the  series  so- 
lutions are  given;  the  R 
independent  transform 
result  (u  - 1+2?)  describes 
the  behavior  of  u in  the 
asymptotic  entry  region. 
Finally,  on  Fig.  10b  the 
experimental  data  of 
reference  [19]  are  also 
included.  For  0.  2>?_>°.  08 
the  entry  and  numerical 
solutions  with  crossflow 


are  in  good  agreement  with  each  other.  The  flow  dependence  on  ? and  therefore 
independence  of  R is  evident.  The  solutions  with  zero  crossflow  are  some- 
what lower  chan  ttA  other  solutions,  but  are  not  significantly  different.  In 
fact,  it  can  be  concluded  that  the  secondary  motion  has  only  a minor  effect  on 
the  axial  behavior.  As  seen  on  Fig.  10b,  agreement  with  the  data  is  good. 

Fully  developed  conditions  are  achieved  for  ? = ?ejJt  s0.  9-1.  0.  In  previous 
studies,  [1-6]  the  entrance  length  has  been  estimated  to  be  in  the  range 
?.nt  = 0.  89-1.  0.  From  Fig.  10b,  it  can  be  inferred  that  the  entry  region  an- 
alysis is  valid  up  to  ? = 0.  2. 

On  Fig-  11  the  axial  velocity  is  compared  with  experimental  data.  Many 
additional  comparisons  were  made,  but  since  the  agreement  is  similar  to  that 
of  Fig.  11,  and  since  the  axial  behavior  has  been  described  quite  well  in  other 
investigations,  these  results  are  not  included  here.  The  axial  velocity  profiles 
along  the  centerline  (z=0)  are  shown  on  Fig.  12. 


( 


Fig-  10a  Centerline  velocity  near  inlet 
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FIG  lot  CENTEALWC  VtLOCITY  VSSUTION 


Fig.  10b  Centerline  velocity  variation 


Fig.  11  Streamwise  velocity  profile 


Profiles  for  the  secondary 
velocity  near  the  corner 
region  are  given  on  Fig.  13. 
Reverse  flow  profiles  have 
been  calculated  in  the 
asymptotic  entry  region. 
Similar  behavior  was  pre- 
dicted by  the  entry  analysis 
of  Part  I,  section  I.  3. 

The  decay  of  the  secondary 
motion  as  fully  developed 
flow  conditions  are  attained 
is  apparent.  The  resolu- 
tion in  the  boundary  layer 
region  is  insufficient  to 
make  quantitative  com- 
parisons of  the  numerical 
and  entry  analysis  for  the 
secondary  motion.  The 
flow  reversal  occurs  over 
one-tenth  of  the  channel 
cross-section  or  one-fifth 
of  the  total  boundary  layer 
thickness  for  ? = 0.2.  For 
the  present  calculations; 
this  is  at  most  two  grid 
points.  Even  with  the 
coarse  mesh,  the  axial 
profiles  (Fig.  14)  agree 
quite  well  with  the  second- 
order  boundary  layer  re- 
sults present  in  Part  I. 
Near  the  end  of  the 
asymptotic  entry  flow 
(Fig.  14b)  large  differences 
near  the  boundary  layer 
edge  are  apparent.  Near 
the  surface  the  agreement 
is' still  quite  good.  In 
Part  I the  extent  of  validity 
of  the  second-order 
boundary  layer  theory  is 
discussed  in  greater  de- 
tail. 

Finally,  the  shear 
parameter  u is  presented 
on  Fig.  15  al^ various  lo- 
cations on  the  surface, 
y»-l.  For  very  small 
values  of  % , the  numerical 
solutions,  the  undisturbed 


Hasius  results  and  the  second-order  boundary  layer  solutions  are  in  close 
agreement.  As  ? increases  the  wall  values  (z=0.95  or  «=  0.90)  decrease 
due  to  corner  layer  interference.  As  the  boundary  layers  grow  further,  the 
mid-point  values  (z  = 0.  5)  also  fall  below  the  centerline  (z  = 0)  solutions. 
Second-order  boundary  layer  theory  is  a significant  improvement  over  the 
Blasius  values  in  the  entry  flow  (?  < 0.  2);  even  for  larger  ? values,  where 
the  theory  Is  questionable,  the  results  are  within  20%  of  the  numerical  so- 
lutions. Once  again,  the  solutions  obtained  by  neglecting  the  secondary 
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Fig.  14b  Streamwise  velocity  Profile:  § = 0.  205,  z=0.  0 


Fig.  15  Shear  stress  parameter:  y:  -1 


motion  are  quite  good.  This 
reinforces  our  earlier  con- 
clusions that  the  axial  flow  is 
virtually  insensitive  to  the 
secondary  motion.  On  the  other 
hand,  the  reverse  is  not  true  as 
the  secondary  motion  near  the 
walls  changes  significantly  with 
mesh  refinement,  which  has  a 
much  smaller  effect  on  the  axial 
flow  behavior.  Accurate  sec- 
ondary flow  solutions  near  the 
walls  would  require  much  finer 
meshes.  Since  there  is  no 
available  data  for  the  secondary 
motion,  finer  mesh  calculations 
were  not  considered.  The  so- 
lutions in  the  potential  core, 
however,  are  reasonably 
accurate.  Isovels  and  Isovors 
are  given  on  Fig.  16.  The  pre- 
dicted core  velocity  behavior 
u - 1+2?  and  v -y  is  closely 
approximated.  The  constant 
vorticity  lines  clearly  show  the 
extent  of  the  boundary  layer 
growth. 


II.  4 Summary 


Numerical  solutions  for  the  flow  in  a square  channel  have  been  obtained 
and  agree  reasonably  well  with  the  results  of  the  entry  region  analysis  of 
Part  I.  The  numerical  procedure  is  also  applicable  for  rectangular  ducts, 
and  provides  solutions  for  the  axial  and  secondary  motion  that  are  fully  con- 
verged at  each  axial  location.  The  calculations  are  Reynolds  number  in- 
dependent, so  that  for  a given  geometry  only  one  set  of  results  is  required 
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for  all  Reynolds  numbers. 

The  secondary  motion  shows 
evidence  of  velocity  reversal 
similar  to  that  predicted  by 
the  entry  analysis.  Although 
the  coarse  mesh  considered 
here  appears  to  be  adequate  for 
the  axial  flow  and  the  second- 
ary motion  in  the  potential 
core,  a finer  mesh  is  re- 
quired to  resolve  the  second- 
ary boundary  layer  behavior 
near  the  walls.  It  is  shown 
that  the  crossflows  have 
little  effect  on  the  axial  flow 
and  that  in  the  entry  region, 
the  latter  is  described 
reasonably  well  with  second- 
order  boundary  layer  theory. 


Fig-  16  Isovels  and  Isovors:  ? = 0.  225 
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